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Abstract 

We study thermodynamic evaporation of Schwarzschild-de Sitter black holes 

in terms of a low energy perturbation theory. A small black hole which is 

far from the cosmological horizon and observers at the spacelike hypersurface 

where black hole attraction and expansion of cosmological horizon balance 

exactly are considered. In the low energy perturbation, scalar field equations 

are solved in both regions of the hypersurface and scattering amplitudes are 

derived. And then the desired thermal temperatures from the two horizons 

are obtained as a "minimal" value of the statistical thermal temperature, and 

the fine-tuning between amplitudes gives a relation of the two temperatures. 
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I. INTRODUCTION 



Recent phenomenological observations show that the expansion of our universe is accel- 
erating [0], which imphes that our universe will look like de Sitter (dS) spacetimes that have 
energies dominated by a negative pressure such as cosmological constant or quintessence. A 
model of black holes corresponding to this universe can be considered, which is described by 
the Einstein-Hilbert action with a positive cosmological constant, and the static and neutral 
solution in this system has been known as a Schwarzschild-de Sitter(SSdS) black hole. In 
this family of black holes, the size of black hole varies from zero to its cosmological limitation 
and then the spacetimes where we live are confined within the cosmological horizon. 

Recently, Hawking et. al. have studied black holes with the cosmological horizon in Refs. 
0-^. If the size of the black hole is much smaller than that of the cosmological horizon, 
the radiation coming from the cosmological horizon is negligible compared to that from the 
black hole horizon, and the effective radiation experienced by observers is similar to the 
case of Schwarzschild(SS) black holes since the thermal temperature from black hole horizon 
Tbh is greater than that from cosmological horizon Tds- So the black hole will radiate more 
than what it receives, and there exists an energy transfer from the black hole horizon r/j to 
the cosmological horizon r^. Finally, the black hole horizon becomes to shrink completely 
while the cosmological horizon grows faster until the black hole disappears to evaporate out. 
Furthermore, the thermodynamic instability for this case has been investigated in Ref. P|. 
For the black hole which grows up to the same size of the cosmological horizon, the radiation 
coming from cosmological horizon equals to the black hole radiation, and it will be in the 
thermal equilibrium. For the nearly degenerate case, if the black hole horizon is shrunk from 
its equilibrium value, the size of black hole increases by "anti-evaporation" while the black 
hole evaporates when we choose a different type of initial perturbations which corresponds 
to a kind of "push" in the radiation bath. 

On the other hand, it would be interesting to study scattering amplitudes of the scalar 
field on the SSdS black hole background since the greybody factor or the decay rate is 
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closely related to the thermal temperature of black holes. Some studies on these quantities 
in terms of the low energy perturbation method have been done for the four- dimensional 
Kerr-Newman(KN) black holes and the Kerr-Newman-de Sitter(KN-dS) black hole 
without considering observers |p. 

In this paper, we shall study on the scattering amplitudes of the massive scalar field on 
the SSdS black hole background by using the low energy perturbation theory. We consider 
a small black hole whose size is much smaller than that of the cosmological horizon and 
assume that the energy of probing fields cu, is negligible compared to the inverse of each 
horizon. In this configuration, the whole spacetimes can be approximately divided by two 
parts - SS black hole and dS spacetimes, and then the massive scalar field equation is 
easily solved on each background. Finally, the scattering problem in SSdS background 
can be solved by matching coefficients between asymptotic and near horizon solutions, and 
defining refiection coefficients. In addition, the temperatures from the two horizons can be 
obtained as "minimal" values of statistical thermal temperatures derived from the reflection 
coefficients. In Sec. II, we present a model of Einstein-Hilbert action with the negative 
cosmological constant, which shows two characteristic properties - SS black hole and dS 
spacetimes, and set up the static configuration of this model. Section III and IV contain 
the specific calculation of low energy perturbation in the attraction dominant region and 
the expansion dominant region. As a result, the expected thermal temperatures at both 
regions are obtained by using the appropriate boundary condition and coefficient matching 
procedure. In Sec. V, some discussions on the fine-tuning of amplitudes between different 
two regions and thermodynamic behaviors of non-degenerate and degenerate cases are briefly 
presented. 

II. MODEL SETUP 

We start with a four-dimensional Einstein-Hilbert action with a cosmological constant 
A > 0, which is given by 
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c(A) _ 

" IQtxG 



^— J d'x^g (R-A), (1) 



where G is a gravitational constant and we set G = 1. From Eq. ([T|), the neutral and static 
solution for the spherically symmetric Einstein equation is described by the SSdS black hole 
metric 

(dsf = -U{r)dt^ + U{r)-^dr^ + r^dVlf^), (2) 

where 

(3) 

r 6 

and is a unit line element of two-sphere, and M is a black hole mass. Note that two 
horizons Vh and Vc {jh < Tc) are only possible for < M < (3V^)~^. In this spacetime, 
there exits a two-spherical hypersurface where the gravitational attraction and cosmological 
expansion balance out exactly, so a geodesic orbit is defined as the position where observers 
need no acceleration in order to stay. From Eq. (^, the geodesic orbit is located at = 
(3MA~^)^/^ with a normalized Killing vector at this hypersurface. Note that it is similar to 
the case of an observer who lives at infinity for the SS black hole. 

We can divide our whole spacetimes into specific four regions as shown in the Fig. (p. 
With the assumption of « Vg « r^, as r goes to the black hole horizon r^, the SS black 
hole metric is prominent in the attraction dominant region (region (I) and (II)) in as 

U(r) ^ 1 - — , (4) 
r 

while if we are in the expansion dominant region (region (III) and (IV)) which is accom- 
phshed by r Tc >> Tg, the metric (^ behaves as hke a dS metric 

f/(r) ^ 1 - |r2. (5) 

Note that the metric element (^) is a patch of Eqs. (§) and (^ at the geodesic orbit Vg 
but this picture is only valid for << r„ << r^. For each region, a field equation can 
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FIG. 1. Configuration of Schwarzschild-de Sitter black holes: Observers are located at the 
two-spherical spacelike hypersurface where cosmological expansion and black hole attraction balance 
out exactly so that there is no acceleration on this hypersurface as like observers at Minkowski 
spacetimes, and as for a thermal radiation, observers will experience thermal temperatures Tbh 
and Tds coming from two horizons. 



be solved and the corresponding amplitudes of solutions between the four regions will be 
matched exactly. Then we define a reflection coefficient as a ratio between ingoing and 
outgoing fluxes in region (II) and region (III), and finally obtain the scattering amplitudes 
and thermal temperatures coming from these two horizons. 

The djTiamical behavior of a massive scalar field under the background (0) is described 

by 

(□ -m^) ^(t,r,^,0) = 0. (6) 

Using the separation of variables, \l/(t, r, 6', 0) = e~*'^*-R(r)B(6')e*^'^, the field equation (^ 
becomes 

{r-U{r)d.R{r)) + - - R{r) = 0, (7) 

J—deismedeQ{e)) - -^^(9) = -u{u+l)Qie), (8) 
sm sm d 

where u is an angular mode. 
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III. LOW ENERGY PERTURBATION IN THE BLACK HOLE REGION 

A. Solution near the black holes(Region (I) and (II)) 

At first, let us consider the field equation (|^) in the near black hole horizon region(Region 
(I) in Fig. (|I])). As r goes to r/i, the metric (0) becomes to be a SS metric, t/(r) ^ (1— 2M/r). 
So the equation of motion ([^ can be written as 

(l - a^Rir) + (? - r|) a.mr) + - ^ - mr) ^ 0. (9) 

where the horizon 2M as r ^ r^. Taking a coordinate transformation such as ^ = 

(r — rh)/r, where z spans from to 1 as r goes from Th to r^, where we assumed that 
Th « Tg, the radial equation becomes 

z{l - z)dlR{z) + (1 - z)d,R{z) + - ^^YTT^ - "^'^') ^(^) = 0- (10) 

Using R{z) = z"^{l — zYg{z) in order to remove two singular points at z = and z = 1, Eq. 
(|T0|) is written as 

z{l - z)dlg{z) + (1 + 2a - (1 + 2a + 2^)z)d,g{z) 

and we have a = ±iujrh and /3 = z/ + 1. Then Eq. (0) becomes 

2(1 - z)dlg{z) + (1 + 2a - (1 + 2a + 2p)z)d,g{z) - ((a + [5f + mV2)^(z) = 0. (12) 

Note that we take the plus signature of a since the solution will be symmetric for changing 
a into —a. The solution of Eq. ([l^) is given as 



R{1) {r) = Cout2'"(l - zfF{a^ + /?, a+ + /5, 1 + 2a; z) 

+QnZ""(l - ^)^F(-a+ + /?, -a_ + /?, 1 - 2a; z), (13) 

where a± = a ± imr/i = i(c<j ± m)rh and Cin and Cout are ingoing and outgoing coefficients, 
respectively. 
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Next, from the solution (pTSD, we shall use a z 1 — z transformation 



{-iyr{a + b-l) ^ (a)„(6), 



E 



1-z) 



T{a-l)T{b-l) fo^Kn + iy. 
X [ln(l - z) - ip{n + 1) - ip{n + / + !) + ipia + n) + %l){b + n)] , (14) 

where ip{z) is a digamma function and a = ±a=p + (3, b = ±a± + (3, and / = 2z/ + 1 in our 



case. As 2 —i> 1, for the lowest order of n, the solution (|l^) is transformed to the region (II) 
by using Eq. ([T^ as 



r(2/3-i)r(i + 2a) , ^ r(2/3-i)r(i-2a) 



+ 



r(a_ + /5)r(a+ + /3) "T(-a+ + /3)r(-a_ + /3) 

P°"*r(l + 2a) ^ Pjfr(l - 2a) 



1-/3 



r / 



'°'''r(l + a_ - P)T{1 + a+-P) ' ^'°r(l - a+ - /3)r(l - a_ - /3) 
where Pjf = -^(1) - ip{2p) + 7/;(-a+ + P) + V^(-a_ + /9) = 1^:°"* and note that Eq. (|15D 



is symmetric for changing (3 into 1 — /?. 



B. Asymptotic solution at Region (II) 

We now consider the asymptotic limit of r — > >> r/j which describes the region (II) 
in Fig. ([^). In this limit, the field equation (^ becomes 

d^,R{r) + ^Rir) + _ "^^^^^^ R{r) = 0, (16) 

where p = ^/cu^ — m^, and the solution is given as Bessel functions, 

^?(r) = ^J_i(,_2,)(pr) + ^J.(,_,,)(pr), (17) 

where Bi and B2 are amplitudes of the solution. The solution (|17D can be expanded by 
keeping the lowest leading order of r and defining new coefficients Sin and i?out as 

Bl = -Bin + -Bout 

52=2(5in-5out), (18) 
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in order to separate wave modes into in-and out-going parts, the solution is written as 



where = (ip) ^^Vr (| + /?) and Q 



1 \ 5(1-2/3) 



Q^^^* = Q^^' and ' is exphcitly calculated as [Q^j^^'Q 



{+)o{-)^ 
13 i 



(19) 

/r (I -/?). Note that 

2 cos /37r _ 2(^1)^ 



7r(l-2/3) 7r(l-2/3)" 



C. Matching process and temperature 



As easily seen from Eqs. (|T5D and (0), by matching these solutions, we obtain the 
following relations. 
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,1-/3 



^ r(2/?-l)r(l-2a) r(2/3 - l)r(l + 2a) 



^(+) ^ -r(-a+ + /3)r(-a_+/5) 



— z 



r(l - 2a)Pjf 



r(l-a+-/3)r(l-a_-/3) 



r(a_ + /3)r(a+ + /3) 

r(i + 20)1); 



lOUt 



out 



5„ 



1-/3 

/3 / 



' r(2/?-l)r(l-2a) ,^ 

'-^in:;^";^ „ . "t L-out 



T(-a;+ + /3)r(-«_ + /5) 
r(l - 2a)Dj: 



r(l + a+ - /3)r(l + a_ - /3) 
r(2/3-l)r(l + 2a) \ 

r(a_ + /3)r(«+ + /?)yl 

r(l + 2a)r'°"* 



(20) 



'°''*r(i + a+ - /?)r(i + a.-l3)^ 

At this stage, we impose a boundary condition Cout = 0, which implies that there are 



no outgoing modes of scalar fields since all modes of fields are absorbed into the black hole 
horizon. Defining a flux as 



= —U{r) {R*{r)drR{r) - R{r)drR*{r)) 
i 



(21) 



we obtain them near the geodesic orbit Vg as 



■out 



-SIR 



out 



|Anr(-l)^r 



/3^-2 



(22) 



The reflection coefficient near the geodesic orbit is defined as a ratio between ingoing and 
outgoing fluxes, and consequently it gives 
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n 



BH 





2 


-Bin 






-Bout 



(23) 



Furthermore, statistical-thermodynamic temperature is defined as a vacuum expectation 



value of a number operator and it is related to the reflection coefficient |]TO|,|TT 

1 n 



< 0\N\0 >-- 



BH 



e'^/Tstat _ 1 1 _ T^gjj ■ 

Therefore, we have a thermal temperature of black holes. 



(24) 



CO 



stat 



InT^BH 

From Eq. (|20|) with the boundary condition Cout 
becomes 



(25) 



0, the logarithm term in Eq. (p5| 



ln7^ 



BH 



In 



-2ii + o{e 



(26) 



where 



i = — r/,(-l)'^+ir(2/5) sin7ra+sin7ra„(r't" -^^;;*") 



27r 



2-2/3 



X 



+ 



,2/3 



in 1 2 



-1 



(27) 



r(l + a+ - f3)T{l + a_ - /?)r(l - «+ - /3)r(l -a^-f3) 

For simplicity, considering the main contribution of s-modes (vanishing angular potential) 
when /5 = 1 as shown in Ref. [f7i|l2|, Eq. (p7|) is simplified as ,^ = —n^iourh/S by using the 
Schwarz's inequality. So, we obtain the "minimal" thermal temperature. 



^min 
stat 



1 



■BH 



1 



(28) 



2'Krh Vtt/ 27rr/i 

where it agrees with the well-known Hawking temperature for the SS black hole. Therefore, 
we have shown that the desired Hawking temperature can be reproduced as a "minimal" 
statistical temperature in attraction dominant region (Region (I) and (H)). 
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IV. LOW ENERGY PERTURBATION IN THE DE-SITTER REGION 

A. Solution near the de Sitter Spacetimes (Region (III) and (IV)) 

Now we study the solution in the expansion dominant region (region (III) and (IV)). As 
r — > Tc, Eq. (|^ becomes 



(rl - r^) df.R{r) + f - - ^ ) drR{r) + 



- R{r) = 0, (29) 



where is nearly (3A ^)^/^ as r goes to Vc- By choosing a coordinate asy = 1 — r^/r^ where 
y spans from to 1, the radial equation (|29| ) is written as 



cu^r^ u^u + 1) m^r^' 



y{l - y)dlR{y) + -(2 - ^y)dyR{y) + 



R{y) = 0. (30) 



2^- ■ y Ay 4(1-1/) 4 ^ 

Singular points at ?/ = and y = 1 can be eliminated by choosing R{y) = y°'{l — y)^u{y), 
and Eq. (|30D is written as 

y{l - y)dlu{y) + (^1 + 2a - ^ (45 + 4/3 + s) y) 9,^i(y) 

m(2/) 



,2^.2- 

c 



1 /3 K^+l) 

2 4 , 



(a + /3)(a + ^ + ^) + ^ 



w(l/) = 0. (31) 



Determining a = iujrc/2 and (3 = v jl = —(1 — /3)/2, Eq. ( PT| ) becomes 
1/(1 - y)dlu{y) + (l + 2a - 1 (4« + 4/3 + s) dyu{y) 



(a + /3)(a + /? + -) + ^ 



By defining 7-1- = (3 ± w9 — 4m^r2)/4, the solution is given as 



uiy) = 0. (32) 



R{y) = AnZ/"(l - yfF{^_ + a + P,^+ + a + ^A + 2a;y) 

+Aout2/~"(l - yfF{^^ -a + (3n+-a + P,l-2a;y), 



(33) 



where 7+ + 7_ = 3/2. Note that the solution (^) follows a z ^ 1 — z transformation given 
by 
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^/ , N r(a)r(c — a — 6) ^ , , , , , X 

F a, b, c; z) = ) > ^ -(F a, 6, a + 6 - c + 1; 1 - z 

r(c — a)T[c — b) 

+ (1 - ^)^-^-^ ^(")^^(^^^+^p) F(c - a, c - 6, c - a - 6 + 1; 1 - z). (34) 

Therefore, for the hmit oi y ^ 1 (r ^ << r^), we obtain the solution in region (III) of 
Eq.(|3D, 



Ry^i{r) 

Anr(i + 25)r (-2/3 - i) ^ Autr(i - 25)r (-2/3 - 1 



r(l + a - /5 - 7_)r(l + a - /3 - 7+) 1(1 - a - /5 - 7_)r(l - a - /5 - 7+) 



r 

Anr(l + 2a)r(2^+i) Ao„tr(l-2a)r(2/3+£ 
^ ' r(a + /3 + 7-)r(6 + /3 + 7+) ^ r(-a + /3 + 7-)r(-a + /3 + 7/ 



B. Asymptotic solution at region (III) 

On the other hand, for the asymptotic hmit r ^ Vg, the equation of motion (^) can be 



written as the same form with Eq. (16) and the solution is given by Bessel functions with 
different amplitudes, Di and D2, 

4?(r) = ^J_.(i_2,)(pr) + ^J.(,„2,)(pO, (36) 



where p = ^uP' — rrP- and /5 = + 1. In the same way for the previous calculation, by 
defining new coefficients D-^^ and -Dout as 

Di = Din + F)out 

I^2 = ^(An-/^out), (37) 



the solution ( PBj ) can be written as 

_ 



4?(r) = An <^-T-^ + ) + ^out <^-TZ7^ - ) • (38) 
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C. Matching process and temperature 



Comparing Eq. (|35[) with Eq. (^) to perform a coefficient match, we obtain the relations 
between coefficients at region (III) and (IV), 



1 



Anrf-^r(l + 25)r(-2/3 



2<) 



A,trf-^r(l-2a)r(-2^-i) 
r(l + 5-/3-7_)r(l + a-/3-7+) ■ r(l-a-/3-7_)r(l-a-/3-7. 



+ 



— I 



Anr-^r(i + 2a)r(2/3 + i; 



+ 



Aoutr-'3r(i-2d)r(2/3 + i; 



2^]|3"^ \ r(a + ^ + 7-)r(d + /3 + 7+) r(-a + /3 + 7_)r(-a + /3 + 7+) 



Anrf-^r(i + 2a)r(-2^ 



21] 



+ i 



(+) 

1 



+ 



A,trf-ir(l - 2(5)r(-2/3 



r(l + a - /3 - 7_)r(l + a - /? - 7+) r(l - a - /? - 7_)r(l - a - /? - 7+) 



Anr-'^r(i + 2«)r(2/3 + i) 



2n(-) 



+ 



Autr-^r(l-2«)r(2/3 + l) 



r(a + + 7-)r(a + /3 + 7+) r(-a + + 7_)r(-5 + /3 + 74 
Now we impose a boundary condition Aout = 0, which means that there are no wave 
modes falling into the cosmological horizon, and we define a reflection coefficient at the 
geodesic orbit T g cLS 



(39) 



n 



DS 







-Dout 


in 







(40) 



where the ingoing and outgoing fluxes at are obtained by the definition (pi]) 

= 8|Anp(-i) V 



= -8|/^outP(-l) V- 



(41) 



With the help of Eq. (p9D, the reflection coefficient (|iO| ) is written as 

1 + 



7^Ds 



(42) 



where 



^ = Tc sin 2a7r sin' 



TT/3 



7^^2(1-/3) (ip) 1-2/3 



X 



+ 



,r(i - a - /3 - 7-)r(i + 5 - /3 - 7-)r(i - 5 - /3 - 7+)r(i + 7^ 



rf7r2(ip)-{i-2/3) 



r(7_ - a + /?)r(7_ + a + /5)r(7+ - a + /?)r(7+ + a + P), 
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(43) 



By using the relations (3 = —(1 — f3)/2 and r(s)r(l — s) = vr/sinvrs, and the thermal 
temperature is given as, 

Tstat = - -^7 (44) 

In 7^DS 

in a leading order of uj. To be explicit, for 1^ = 0, i.e. /? = 0, Eq. ( ^31) is simplified as 

i= , , , , . (45) 

y COS"' 77(7+ — 7_) — sin 27ra 

Expanding Eq. ( ^5|) in a low energy limit, ^ = 27ra; + 0{m? ^ a^) and the "minimal" thermal 
temperature can be obtained as 

^2? = ^Ds ^ (46) 
which agrees with the thermal temperature of the dS spacetimes. 



V. DISCUSSIONS 

In summary, we have studied the scattering problem of the massive scalar field on the 
SSdS background by assuming the low energy of the probing field. For the small static SSdS 
black hole, the whole spacetimes can be considered as a patch of two different regions (SS 
black hole and dS spacetimes) and the geodesic orbit where the black hole attraction and 
cosmic expansion balance out is introduced. Using the low energy approximation, the corre- 
sponding statistical temperatures from two horizons are derived, and it is interesting to note 
that the minimal bounds of these temperatures are equivalent to the thermal temperatures 
obtained from the metric. 

As a comment, we should perform a fine-tuning process between the wave functions at 
the region (II) and (III) in Sec. II. The wave functions of the massive scalar field should 
be connected smoothly at the geodesic orbit since it is regarded as a continuous function 
in the whole regions. Therefore, from Eqs. (|T9|) and (|38D, we define i?in = Ni^Di^ and 
Bout = A'out-Dout, respectively, where A^in and N^nt are proportional constants. And then, 
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Eqs. (H) and (ID produce a relation i?DS^BH = by defining |Qp = \N-J'^ /\Nout\'^ 

and it determines, 

= e'^(^sH+'^Ds)^ (47) 

which is a fine-tuning for mismatched wave functions at the geodesic orbit and it correlates 
the cosmological horizon Tc with the black hole horizon r/j. By means of this final coefficient- 
tuning, the wave functions of probing fields are continuous and connected smoothly in the 
whole area. 

Now let us discuss about the qualitative behaviors of non-degenerate and nearly degen- 
erate cases of SSdS black holes. As mentioned above, our whole calculations have been done 
by assuming that the black hole horizon size is much smaller than that of the cosmologi- 
cal horizon, which describes a small SSdS black hole. Since the black hole will be hotter 
than the cosmological horizon, the one will evaporate faster than the other, and the black 
hole horizon will disappear at the end of evaporation. And then the SSdS spacetimes may 
become the pure dS spacetimes at the end of the evaporation. At this stage, one may ask 
a question : "After evaporation at final stage, is it really possible to transit from the SSdS 
spacetimes to the dS ones?" In fact, to answer this question explicitly, one should study this 
model with the gravitational back reaction by using the time evolution of the small SSdS 
black hole, however we do not discuss anymore because of complexity of this analysis. 

On the other hand, we naturally consider the degenerate case of horizons that the black 
hole horizon is equal to the cosmological horizon. In this case, we expect that the black hole 
will be in the thermal equilibrium state since the black hole and cosmological horizons have 
the same surface gravity. Therefore, when the black hole horizon (r/^) approaches to the 
cosmological horizon (rc), the two 5*^ hypersufaces (r/j and Tc) meet each other in a hyper- 
surface (r^), and this is somewhat trivial to discuss. To investigate the small perturbation 
from the thermal equilibrium, we can consider the nearly degenerate case whose metric is 
given by 

{dsf = ^ (l + -ScosR) sin^Rd'^T + ^ (l - -ScosR) d^R + ^ (1 - 25cos) ^^^(2) (48) 

AV3 / A\3 / A 
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with defining new time and radial coordinates T and R by 

1 „ 

t 



^ ^1 - 5cosR - , (49) 



5VA 
L( 

VaV" 6 

and 9M^A = 1 — 35^, where < 5 < 1. And the temperatures from two horizons are given 
by Ref. i 

rh,c = ^v^(l±^5)+0(52), (50) 

where the upper (lower) sign is for the black hole (cosmological) horizon. Note that if 5 goes 
to zero, Eq. ([50| ) describes the degenerate temperature characterized by only a cosmological 
constant A, and the instability of the nearly degenerate SSdS black hole has been studied 
by adding the one loop quantum correction terms As a further work, if the scattering 
problem in the nearly degenerate SSdS black hole background is solved by use of the low 
energy perturbation method, the results of both temperatures will be compatible with Eq. 



Acknowledgments 

This work has been achieved partially for participating in the APCTP-PIMS Summer Work- 
shop at SFU, Canada. W. T. Kim and J. J. Oh would like to thank PIMS for hospitality 
during staying at SFU. J. J. Oh and K. H. Yee would like to thank V. P. Frolov for helpful 
comments while they were participating in the 6th APCTP Winter School at POSTECH, 
Korea. This work was supported by the Korea Research Foundation Grant, KRF-2001-015- 
DP0083. 



15 



REFERENCES 



[1] S. Perlmutter et al, Measurements of the Cosmological Parameters Q and A from the 
First 7 Supernovae at z > 0.35, Astrophys. J. 483, 565 (1997) , [ Jastro-ph/9608192| ] ; 



S. Perlmutter, Supernovae, dark energy, and the accelerating universe: The status of 
the cosmological parameters, in Proc. of the 19th Int. Sym. on Photon and Lepton 
Interactions at High Energy LP99, Int. J. Mod. Phys. A15S1, 715 (2000). 

[2] R. Bousso and S. W. Hawking, Pair creation of black holes during inflation, Phys. Rev. 
D54, 6312 (1996). 

[3] R. Bousso and S. W. Hawking, (Anti-) evaporation of Schwarzschild-de Sitter black holes, 
Phys. Rev. D57, 2436 (1998). 

[4] G. W. Gibbons and S. W. Hawking, Action Integrals and Partition Functions in Quan- 
tum Gravity, Phys. Rev. D 15, 2752 (1977). 

[5] P. Ginsparg and M. J. Perry, Semiclassical Perdurance of De Sitter Space, Nucl. Phys. 
B222, 245 (1983). 

[6] J. Maldacena and A. Strominger, Universal low-energy dynamics for rotating black holes, 
Phys. Rev. D56, 4975 (1997). 

[7] M. Cvetic and F. Larsen, Greybody factors for black holes in four dimensions : Particles 
with spin, Phys. Rev. D57, 6297 (1998). 

[8] H. Suzuki, E. Takasugi, and H. Umetsu, Absorption rate of the Kerr-de Sitter black 
hole and the Kerr-N ewman-de Sitter black hole. Prog. Theor. Phys. 103, 723 (2000) ; 
Analytic Solutions of Teukolsky Equation in Kerr-de Sitter and Kerr-Newman-de Sitter 
Geometries, Prog. Theor. Phys. 102, 253 (1999). 

[9] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, Dover Publi- 
cation Inc., New York, ninth printing, (1970). 



16 



[10] W. T. Kim and J. J. Oh, Dilaton driven Hawking radiation in AdS2 black hole, Phys. 
Lett. B461, 189 (1999). 

[11] K. Ghoroku and A. L. Larsen, Hawking temperature and string scattering off the 2+1 
black hole, Phys. Lett. B328, 28 (1994). 

[12] V. P. Prolov, in private communication for participating in "6th APCTP Winter School 
on Black Hole Astrophysics 2002" at POSTECH, Korea. 



17 



